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Pauli-Villars rcgularization is successfully applied to nonperturbative calculations 
in a (3 + l)-dimcnsional light-cone model. Numerical results obtained with dis- 
cretized light-cone quantization compare favorably with the analytic solution. 



1 Introduction 

One of the challenges of using light-cone methods 0! to solve nonperturbative 
problems in field theories such as quantum chromodynamics (QCD) is to de- 
velop a nonperturbativ^Efigularization and renormalization scheme. Here we 
discuss a new approach&Q based on generalizecLPauli-Villars regularizationH 
and discretized light-cone quantization (DLCQ)l£l The DLCQ Fock basis is ex- 
panded to include momentum states of Pauli-Villars particles which provide 
the cancelations needed to regulate ultraviolet divergences. The determina- 
tion of the number of Pauli-Villars particles and of their coupling strengths 
is done in perturbation theory, where loop integrals are rendered absolutely 
convergent. 

We test this approach on a simple model designed to have an analytic 
solution for the lowest massive state. The model is related to an equal-time 
model of Greenberg and Sphweber0 which has also been translated to the light 
cone by Glazek and PerryO 

2 Light-Cone Coordinates 

We define light-cone coordinates and momentum components by0 

x =t±z, x± = (x,y), (1) 



a To appear in the proceedings of the third workshop on Continuous Advances in QCD, 
Minneapolis, Minnesota, April 16-19, 1998. 



and 

p ± = E± Pz , PX = (jPx,Py), P=(p + ,P±_)- (2) 

The dot product is written 

P ■ x = ^{p + x~ + p~x + ) - p± ■ xj_ . (3) 

The time variable is taken to be x + . The conjugate variable p~ is the light- 
cone energy, and the associated operator V~ determines the time evolution 
of the system. The mass-squared operator, frequently called the light-cone 
Hamiltonian0 is 

#lc = V+V- - V\ , (4) 

where V + and V± are momentum operators conjugate to x~ and xj_- In these 
coordinates, the fundamental eigenvalue problem is 

H LC ^ = MH, (5) 

where M is the mass of the state "J. nn 

The use of light-cone coordinates has several advantagesEra The genera- 
tors of the Poincare algebra have the largest possible nondynamical subset; in 
particular, boosts are kinematical. The perturbative vacuum is the physical 
vacuum, because p + = \/p 2 + m 2 + p z > 0, and there is then no need to 
compute the vacuum statej^] Fock-state expansions are well defined with no 
disconnected vacuum pieces. 

For field theories quantized in light-cone coordinates one does need to be 
careful about the use of Pauli-Villars regulators. The number of Pauli-Villars 
particles used in ordinary Feynman perturbation theory may not suffice. One 
concrete example can be found in Yukawa theorya To regulate the one-loop self- 
energy three heavy bosons are required, rather than the usual one. The reason 
that one is not sufficient is that symmetric integration techniques cannot be 
employed; the use of such techniques in equal-time quantization amounts to a 
regularization prescription which supplements the Pauli-Villars regularization. 

To see the need for three heavy bosons explicitly, consider the self-energy 
integral 



dl+d 2 l ± 
l+{q+ -l+) 2 



{q+) 2 \\ + (2g+ - l+) 2 M 
M 2 -D 



2 



■0(A 2 -D), (6) 



'For theories with complicated vacuum structure in equal time fnBimulations this simplicity 
of the vacuum is offset by a more complicated operator structurecl 
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where /i is the ordinary boson mass, M is the fermion mass, and D is the 
invariant mass of the intermediate state 



D = 



l+/q+ 



M 2 + \\ 



(7) 



The integral has been made finite by an invariant-mass cutoff A 2 , and can be 
done exactly; however, an expansion in powers of the fermion mass is much 
more instructive. We find 



9 5^ 2 



— /i 2 liiA 2 +/z 2 m/z 2 - ^ 



+ M 31nA z - 31n/z - - + -fV 

2 A 2 



(8) 



The leading term violates the chiral symmetry of the original theory; to remove 
this term from the infinite-cutoff limit requires three Pauli-Villars bosons with 
different masses /itj. The subtracted integral is 



(9) 



i=l 



and the Ci are chosen to satisfy 

3 3 



l + ]Ta = 0, M 2 +^Q M 2 = 0, ^a M 2 ln( M 2 / M 2 )=0. (10) 



i=l 



Given this perturbative analysis of the regularization, the conjecture is that 
fermion self-energies in a nonperturbative calculation will be regularized by 
the same Pauli-Villars bosons with the same coupling strengths Ci|] 



3 A Soluble Model 

To test this idea, we have constructed a model remotely related to Yukawa 
theory but which has the advantage of being analytically soluble. This type of 
model was investigated in equal-time quantization by Greenberg and Schwc- 
berQ and on the light cone by Glazek and PerryQ It can be obtained from the 

c To completely regularize Yukawa theory will require a <f> 4 term and perhaps a Pauli— Villars 
fermion. 
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Yukawa Hamiltonian0 by a number of severe modifications. The momentum 
dependence in the fermion kinetic energy becomes (M§ + MqP + ) / P + , where 
p + and P + are the longitudinal momenta of the fermion and system, respec- 
tively^ Only the no-flip three-point vertex is kept and then only in a modified 
form where the longitudinal momentum dependence is simplified. The struc- 
ture of the chosen interactions results in states populated by a fixed number 
of fermions and a cloud of bosons. We study only the lowest state with one 
fermion. One Pauli-Villars field is found to be sufficient in this case. 
The resulting light-cone Hamiltonian iJ£^ = P + P~ s is given by 

dp + d 2 pj_ 



(7 

dq + d 2 q_ 



-P 1 



q + £°a + q + a k a H 



167r 3 g 

f dp^d 2 p ±1 f dp^d 2 p± 2 f dq + d 2 qx \^ , t , 



(11) 



lQir^Pi 

+ \ 1 



167r 3 j>2 



(§0 4%-h 2 +2) + (|) 



where 



= 16tt 



(12) 



The state vector is written as a Fock-state expansion 



drj~d 2 r±j 



n ni 

<s(p - p - E % - E ^)0 (n,Bl) (2*. £j ; e) 

V?T-!ni! - -i A r l -j 



(13) 



d This t crrn has a structure similar to that of the self-induced inertia term shown in Eq. (C.2) 
of Ref . [lo| 



with n the number of ordinary bosons and n\ the number of Pauli-Villars 
bosons. The normalization of the state is chosen to be 



< f • $> a = 167r 3 P+(5(P' - P) . 
For the boson amplitudes this implies 

^En/ ^Axiii/ dr t d2r ^ 

■ ■ i j 

A solution to the mass eigenvalue problem 



(14) 



(15) 



(16) 



must satisfy the following coupled set of integral equations: 

^Kg^p) 



M 2 Ml M' p+ V 'J—!U Y.' 1 ' ' 1 



+ \ 7 



■\/l67r 3 <7 



(17) 



+ ivni + 1 



dr + d 2 r±_ ( ' p + — r + 



^ ni+1) {g v r_ v r_,p-r) 



i 1 / p + \ 7 

+ V^?^/l1^U + +^J 

'' ni ~ 1) (£ i .Ei.-".Ej-i.Ej+i.---.rni»H + ^)} • 



X 

The solution is 



n 



V / 16^?(M? + ^ J )' 



provided that M — M and 



16tt 2 7 + 1/2 



(18) 



(19) 



Given the structure of the individual amplitudes, a natural choice for the value 
of the interaction parameter 7 is seen to be 1/2, because each amplitude is then 
proportional to the square root of the product of all longitudinal momenta. 
With this value of 1/2, the wave function normalization can be reduced to a 
rapidly converging infinite sum 



00 



1 



(2/M) 2n (.9/Mi) 



2ni 



(2n + 2m + l)!n!m! (167r 2 )"+"i 



(20) 



To fix the coupling g we set the value of the expectation value (:</> 2 (0):) 
:</> 2 (0): <& a . For the analytic solution (with 7 = 1/2) it reduces to 



<:0 2 (O):> = E 



2Zn 



(ff/M) 2n (s/Mi) 2ni 



(2n + 2m)!n!ni! (16n 2 ) n+n ^ 



(21) 



From a numerical solution the quantity can be computed fairly efficiently in a 
sum similar to the normalization sum 



(^(O):) = E II / d tfd 2 q^i[ / dr+d 2 r ±j 



(22) 



n— l,ni= 



^(sw^-Ea-E^) 



with the integrals computed from discrete approximations. 

As a prediction of the model, we compute a distribution function for the 
physical bosons 

n ,. n\ „ n 

Mf) = EIl/ ^Axill/ dr+<Pr ±j '£5{y-qt/P + ) 

r, A J A J A 1 



^(sw^-Ea-E^) 



(23) 
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and the average multiplicity 

(«b) = / fB(y)dy. 
Jo 

For the analytic solution we obtain 

Zny(l - ( fl //i) 2 »( 5 / M i) 2ni 



My) = E 



(2n + 2ri! - l)!n!ni! (167r 2 )"+™i 



and 



\ n B) 



E 



(5/M) 2 "(5/Mi) 



2iii 



(2n + 2ni + l)!n!ni! (167r 2 ) rl +"i 



(24) 



(25) 



(26) 



Another prediction is the slops of the form factor for the dressed fermion. 
An expression can be constructed u from the formalism developed by Brodsky 
and Drellcil For the analytic solution of the model we obtain 



no) = -E 



Zjn/f + m/rt) (g//i) 2 "(g /Mi) 2 " 1 
(2n + 2n x +3)!n!n x ! (16tt 2 ) 



2 \n-\-7li 



Numerically, one can estimate F'(0) from 

n „ Hi n 



r ±3 



(27) 



(28) 



E 



E 



-V^'-^^r^P-Ei-E^) 



which differs from F'(0) by surface terms which vanish as A — > oo. 
4 Numerical methods and results 

We solve the model numerically by using DLCGji We impose periodic boundary 
conditions for bosons and antiperiodic conditions for fermions in a light-cone 
box — L < x~ < L, —Lj_ < x,y < L±. This introduces a discrete grid: 



+ 7T 7T 7T . 

f 1 ^t", P± -> (y— «x, y— %J 



(29) 
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Integrals are replaced by discrete sums 



2 N 

J dp+ J dV/(p + ,Pi) (tz) ^ ^ ! >!•< I L. 11 TT/'Z. _ 



n n x .n v =—N± 



(30) 

TheJamit L — > oo can be exchanged for a limit in terms of the integer resolu- 
tions 

K=-P+. (31) 

Longitudinal momentum fractions are given by x = p + /P + — > n/K, with n 
odd for fermions and even for bosons. The light-cone Hamiltonian -ZZlc is 
independent of L, but does depend on K. 

Because all n are positive, DLCQ automatically limits the number of par- 
ticles to no more than ~ K/2. The integers n x and n y range between limits 
associated with some maximum integer N± fixed by the invariant-mass cutoff 

m i +p ±i < A 2 (32) 



imposed for each constituent. The eigenvalue problem ( |l7| ) is then converted 
to a finite matrix problem. Typical basis sizes are given in Table 0. The Hamil- 
tonian matrix is quite sparse. TheJnwest eigenvalue of the matrix is extracted 
with use of the Lanczos algorithmE3 for complex symmetric matrices.^ 

Some of the results obtained are listed in Table || and displayed in Figs. [I] 
and ||. Figure [j] shows that the results are quite insensitive to numerical reso- 
lution, while Fig. || illustrates how the analytic solution is approached as the 
cutoff is increased. These results demonstrate that the DLCQ approximation 
yields a good representation of the solution to the model. 

5 Summary 

We have developed and used a simple soluble model to test the feasibility of 
Pauli-Villars regularization in DLCQ. The number of Pauli-Villars Fock states 
required is not prohibitive and good results are obtained. From this success 
we can work toward approximation of Yukawa theory itself by adding fermion 
dynamics and gradually reinstating the full complexity of the interactions. We 
can also consider application to other theories, perhaps including QCD. 



e The imaginary couplings of the Pauli-Villars particles make the Hamiltonian matrix com- 
plex symmetric. 
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Table 1: Basis sizes for DLCQ calculations in the soluble model with parameters M 2 = fj, 2 , 
fi\ = 10^t 2 , and A 2 = 50/j 2 . The numbers of physical states are in parentheses. 



N ± 


K 


7 


9 


11 


13 


15 


17 


1 


18 


38 


36 


65 


110 


185 




(7) 


(12) 


(19) 


(30) 


(45) 


(67) 


2 


218 


265 


590 


1120 


822 


1410 




(127) 


(119) 


(343) 


(754) 


(453) 


(626) 


3 


958 


1408 


4460 


17031 


22486 


21635 




(367) 


(736) 


(2671) 


(9230) 


(13213) 


(13531) 


4 


3714 


9259 


49394 


50966 


110254 


328966 




(1399) 


(5913) 


(32363) 


(32124) 


(55319) 


(172247) 


5 


13702 


54100 


95176 


386140 


1553576 






(5699) 


(28065) 


(66371) 


(232400) 


(1038070) 




6 


35666 


126748 


536758 


2907158 








(12991) 


(69245) 


(391511) 


(2107688) 






7 


79794 


519325 


1317392 










(32891) 


(276299) 


(1008539) 








8 


172118 


1165832 












(61947) 


(687394) 











Table 2: Numerical parameter values and results from solving the model eigenvalue problem. 
The physical parameter values were M 2 = for the fermion mass, fi\ = Wfi 2 for the Pauli— 
Villars mass, and (:<fi 2 (0):) = 1 to fix the coupling g. 



(A/m) 2 


K 


N±_ 


LiLj_/ir 


(M //i) 2 


9/» 


(n B ) 


100/^(0) 


50 


11 


4 


0.8165 


0.8547 


13.293 





177 


-0.751 


50 


13 


4 


0.8165 


0.8518 


13.230 





172 


-1.015 


50 


15 


4 


0.8165 


0.8408 


13.556 





178 


-0.715 


50 


17 


4 


0.8165 


0.8289 


13.392 





180 


-0.565 


50 


9 


5 


1.2062 


0.8601 


14.023 





179 


-0.547 


50 


9 


6 


1.2247 


0.8377 


14.323 





179 


-0.582 


50 


9 


7 


1.4289 


0.8302 


14.386 





179 


-0.658 


50 


9 


5 


1.2062 


0.8601 


14.023 





179 


-0.547 


100 


9 


5 


0.7143 


1.0520 


12.565 





174 


-0.239 


200 


9 


5 


0.5025 


1.1980 


10.191 





172 


-0.139 


oo 




ana 


lytic 


1.0000 


13.148 





160 


-0.786 



9 




10 



0.30 




Figure 2: The boson distribution function fs for different cutoff values, with (:</> 2 (0):} = 1 
and numerical resolution set at K = 9 and N± = 5. The solid line is the analytic result at 
A 2 = oo. 
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